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We present a systematic study of the role of Anderson orthogonality for the dynamics after a 
quantum quench in quantum impurity models, using the numerical renormalization group. As shown 
by Anderson in 1967, the scattering phase shifts of the single-particle wave functions constituting the 
Fermi sea have to adjust in response to the sudden change in the local parameters of the Hamiltonian, 
causing the initial and final ground states to be orthogonal. This so-called Anderson orthogonality 
catastrophe also influences dynamical properties, such as spectral functions. Their low-frequency 
behaviour shows nontrivial power laws, with exponents that can be understood using a generalization 
of simple arguments introduced by Hopfield and others for the X-ray edge singularity problem. The 
goal of this work is to formulate these generalized rules, as well as to numerically illustrate them 
for quantum quenches in impurity models involving local interactions. As a simple yet instructive 
example, we use the interacting resonant level model as testing ground for our generalized Hopfield 
rule. We then analyse a model exhibiting population switching between two dot levels as a function 
of gate voltage, probed by a local Coulomb interaction with an additional lead serving as charge 
sensor. We confirm a recent prediction that charge sensing can induce a quantum phase transition 
for this system, causing the population switch to become abrupt. We elucidate the role of Anderson 
orthogonality for this effect by explicitly calculating the relevant orthogonality exponents. 

PACS numbers: 02.70.-c, 05.10.Cc, 71.27.+a, 72.10.Fk, 73.21.La, 75.20.Hr, 78.20.Bh 



I. INTRODUCTION 

The Anderson orthogonality (AO) catastrophe^ refers 
to the response of a Fermi sea to a change in a local scat- 
tering potential, described, say, by a change in Hamil- 
tonian from Hi to Hf. Such a change induces changes 
in the scattering phase shifts of all single-particle wave 
fimctions. This causes the initial ground state |Gi) of Hi 
and the final ground state \G[) of Hi, both describing a 
filled Fermi sea but w.r.t. different single-particle wave 
functions, to be orthogonal in the thermodynamic limit, 
even if the changes in the single-particle wave functions 
are minute. The overlap of the respective ground states 
scales as^'^ 

|(Gi|Gf)|^iV-5^Ao, (1) 

where N is the number of particles in the system, and the 
exponent Aao characterizes the degree of orthogonality. 

AO underlies the physics of numerous dynamical 
phenomena such as the Fermi edge singularity,^"^ the 
Altshuler-Aronov zero bias anomaly^ in disordered con- 
ductors, tunnelling into strongly interacting Luttinger 
liquids, ^"^^ slow relaxation in electron glasses, ^^'^'^ and 
optical absorption involving a Kondo exciton,^^"^^ where 
photon absorption induces a local quantum quench, to 
name but a few. Recently, AO has also been evoked^^'^^ 
in an analysis of population switching (PS) in quantum 
dots (the fact that the population of individual levels of a 
quantum dot may vary non-monotonically with the gate 
voltage), and was argued to lead, under certain condi- 
tions involving a local Coulomb interaction with a nearby 



charge sensor, to a quantum phase transition. 

One of the goals of the present work is to analyse the 
latter prediction in quantitative detail. Another is to 
generalize arguments that were given in Refs. 14-16, for 
the role of AO for spectral functions of the excitonic An- 
derson model, to related models with a similar struc- 
ture. Thus, we present a systematic study of the role of 
Anderson orthogonality for the dynamics after a quan- 
tum quench in quantum impurity models involving local 
interactions, using the numerical renormalization group 
(NRG).^^'^° We thereby extend a recent study,^-^ which 
showed how Aao can be calculated very accurately (with 
errors below 1%) by using NRG to directly evaluate over- 
laps such as (Gi|Gt), to the domain of dynamical quan- 
tities. 

The spectral functions that characterize a local 
quantum quench typically show power-law behaviour, 
~ u}~^~^'^^, in the limit of small frequencies, where ry typ- 
ically depends on Aao-'^"^ For the case of the X-ray edge 
singularity, Hopfield* gave a simple argument to explain 
the relation between Aao and rj. We frame Hopfield's 
argument in a more general setting and numerically il- 
lustrate the validity of the resulting generalized Hop- 
field rule (Eq. (25) below) for several nontrivial models. 
In particular, we also analyse how this power-law be- 
haviour is modified at low frequencies when one adds to 
the Hamiltonian an extra tunnelling term, that describes 
transitions between the Hilbert spaces characterizing the 
"initial" and "final" configurations. This effect plays a 
crucial role in understanding the abovcmcntioned quan- 
tum phase transition for population switching. 

The paper is organized as follows. In Sec. II we review 
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various consequences of AO in different but related set- 
tings, and formulate the abovementioned generalization 
of Hopfield's rule. In Sec. Ill we illustrate this rule for the 
spinless interacting resonant level model (IRLM), involv- 
ing a single localized level interacting with the Fermi sea 
of a single lead. We consider this model without and with 
tunnelling, and study a quantum quench of the energy 
of its local level, focussing on signatures of AO in each 
case. Finally, in Sec. IV and Sec. V we discuss population 
switching without and with a charge sensor, respectively, 
confirming that if the sensor is sufficiently strongly cou- 
pled, AO indeed does cause population switching to be- 
come a sharp quantum phase transition. Section VI offers 
concluding remarks and outlines prospective applications 
of the present analysis. 



II. VARIOUS CONSEQUENCES OF 
ANDERSON ORTHOGONALITY 

In this section we review various consequences of AO, 
in different but related settings. We begin by recalling 
two well-known facts: first, the relation between the ex- 
ponent Aao and the charge that is displaced due to the 
quantum quench, Ach; and second, the role of Aao in 
determining the asymptotic long-time power-law decay 
of correlation functions Gx{t) involving an operator 
that connects the initial and final ground state. 

Then we consider the spectral function ^x(w) asso- 
ciated with Gxit): which correspondingly shows asymp- 
totic power-law behaviour, ^ for small frequen- 
cies, where the exponent rj depends on Aao- We recall 
and generalize an argument due to Hopfield, that extends 
the relation between rj and Aao to composite local op- 
erators. Finally, we recapitulate how all these quantities 
can be calculated using NRG. 

For simplicity, we assume in most of this section that 
the Fermi sea consists only of a single species of (spinless) 
electrons. The generalization to several channels needed 
in subsequent sections (in particular for discussing PS), 
is straightforward and will be introduced later as needed. 

Although the concepts summarized in subsections II B 
to II E below apply quite generically to a wide range of 
impurity models, for definiteness we will illustrate them 
by referring to a particularly simple example, to be called 
the "local charge model" (LCM), which we define next. 




Figure 1: (a) Cartoon of the Hamiltonian (2) for the LCM. 
(b) to (g) Cartoons of the occupation of the dot and a half- 
filled lead, for U > 0, for several states discussed in the text, 
(b) and (c) give two equivalent depictions of the ground state 
|Go) of Hq. (c) depicts the fact that \go) can be written as 
a superposition of the form 10)clQ)rcst + !l)c|Q — l)rcst, in- 
dicating complementary occupations of the first site and the 
rest of a half-filled Wilson chain (defined in Sec. TIE below). 
Here |0)c (which obeys c|0)c = 0) and |l)c = c^|0)c describe 
the first site of the Wilson chain being empty or filled, re- 
spectively; the charge in the rest of the Wilson chain is corre- 
spondingly distributed in such a way that both components of 
the superposition have the same total charge, Q. (d) depicts 
the ground state Id) of Hi, indicating that charge on the dot 
pushes charge in the lead away from the dot site, (e) shows 
the effect of applying to |Go), the latter depicted according 
to (b). Similarly, (f) and (g) show the efi'ect of applying c'd^ 
or cct' to I Go), the latter depicted according to (c). The dis- 
placed charge flowing inwards from inflnity towards the dot 
as each of the states (e) to (g) evolves to the final ground 
state |Gi) of (d) is Ad < 0, - 1 < or 4- 1 > 0, re- 
spectively. Comparison of (f) and (g) with (e) shows average 
charge differences of -1-1 and —1, respectively, in accord with 
the Hopfield-type argument summarized by Eq. (15). 



A. Local charge model 

The LCM describes a single spinless localized level, 
to be called dot level (alluding to a localized level in 
a quantum dot), interacting with a single Fermi sea of 
spinless electrons [see Fig. 1(a)]: 

HhCMind) ^Undc^c + ^eclc^ . (2) 



Here and d are annihilation operators for Fermi sea 
states and the dot state, respectively, fid = d^d counts the 
number of dot electrons, and c = ■0(0) = X^e '^e destroys 
a Fermi sea electron at the position of the dot. The 
interaction is taken to be repulsive, U > 0. There is no 
tunnelling between dot and sea. Therefore, the Hilbert 
space separates into two distinct sectors, in which the 
local charge operator hd has eigenvalues Ud — and — 
1, respectively. The Hamiltonians describing the Fermi 
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sea in the two distinct sectors are 

Hq = Hi^cuind = 0) = ^£c)^c^ , (3a) 

HiiU) = 7?LCM("d = l)=5l£4c. + C/ctc. (3b) 
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We will denote their respective ground states [illustrated 
in Figs. l(b,c) and 1(d), respectively] by 

|Go) = |0)|go), |Gi) = |l)|.gi), (4) 

where |0) and |1) = d}\0) describe the dot state with 
charge or 1, respectively, and jf^o) and jf^i) the corre- 
sponding Fermi sea ground states. 

The LCM contains all ingredients needed for AO, hence 
we will repeatedly refer to it below as an explicit example 
of the general arguments to be presented. [Corresponding 
LCM passages will sometimes appear in square brackets, 
so as not to disrupt the general flow of the discussion.] 
Explicit numerical results for the LCM will be presented 
in Sec. Ill A below. 



B. AO and the displaced charge 

For the ensuing discussions, it will be useful to distin- 
guish between two types of quenches, to be called type 1 
and 2, which we now discuss in turn. 

Type 1 quench: For a type 1 quench, some parameter 
of the Hamiltonian is changed abruptly (e.g. by a sudden 
change of gate voltage for one of the gates defining a 
quantum dot). Taking the LCM as an example, suppose 
that the value of the interaction in the LCM is changed 
suddenly from U to U' for a fixed local charge of = 1 . 
This corresponds to a type 1 quench with 

Hi^HiiU), H[^Hi{U'), (5a) 

|Gi) = |l)|5i,i), |Gf) = |l)|gi,f). (5b) 

The overlap of initial and final ground states, 

l(Gi|Gf)| = |(gi,i|gi,f)|^A^-^^Ao, (g) 

will vanish in the thermodynamic limit due to AO, since 
the two Fermi sea states \gi,i) and [51, f) feel scattering 
potentials of different strengths. 

In his classic 1967 paper, Anderson showed that for this 
type of situation the exponent Aao in Eq. (6) is equal to 
the change in scattering phase shifts at the Fermi surface 
divided by tt, in reaction to the change in the strength 
of the scattering potential. According to the Friedel sum 
rule,^^"^^ the change in phase shifts divided by tt, in turn, 
is equal to the displaced charge Ach (in units of e) that 
flows inward from infinity into a large but finite volume 
(say Viargc) surrounding the scattering site, in reaction to 
the change in scattering potential, so that Aao — Ach- 
To be explicit, 

Ach = (GflntotjGf) - (GilntotjGi) , (7) 



where ntot = nsea + 'T-dot counts the total number of elec- 
trons within Viaigc, with hsoa counting the Fermi sea elec- 
trons and ndot counting the electrons on the dot. [For the 
LCM, ndot = fid-] 

The relative sign between Aao and Ach (+ not — ) is 
a matter of convention, which does not affect the orthog- 
onality exponent A^q. Our convention, which agrees 
with standard usage, is such that Aao > (or < 0) 
if the change in local potential induces electrons to flow 
inward toward (outward away from) the scattering site. 

For the LCM quench of Eq. (5) above, the initial and 
final states have the same dot charge, = 1, hence 
the displaced charge reduces to Ach = (<?i,f |^sca|5i,f) — 
(<?i,i|^sca|5'i.i) • However, such a simplification will not 
occur for more complex impurity models involving tun- 
neling between dot and lead [of the form {d^c + c^d)], so 
that the local charge is not conserved. Examples are the 
interacting resonant level model [Eq. (36) below], or the 
single-impurity Anderson model [Eq. (47) below]. 

For such a model, consider a type 1 quench from Hi 
to iJf , implemented by a sudden change in one or several 
model parameters, in analogy to Eq. (5a). Although the 
corresponding ground states |Gi) and |Gt) wiU no longer 
have the simple factorized form of Eq. (5b), they will still 
exhibit AO as in Eq. (1). Moreover, the decay exponent is 
still equal to the displaced charge, Aao = Ach, given by 
Eq. (7). (For a NRG verification of this fact, see Ref. 21.) 

Type 2 quench: For a type 2 quench, all model param- 
eters are kept constant, but the system is switched sud- 
denly between two dynamically disconnected sectors of 
Hilbert space characterized by different conserved quan- 
tum numbers. Taking again the LCM of Eq. (2) as an ex- 
ample, suppose that the local charge is suddenly changed, 
say from rt^ = to 1, while all model parameters are kept 
constant. This corresponds to a type 2 quench with 

ill = Ho Hi = Hi, (8a) 

|Gi) = |0)|5o), |Gf) = |l)|.gi). (8b) 

A physical example of such a quench would be core level 
X-ray photoemission spectroscopy (XPS), where an inci- 
dent X-ray photon is absorbed by an atom in a crystal, 
accompanied by the ejection of a core electron from the 
material. This amounts to the sudden creation of a core 
hole, which subsequently interacts with the Fermi sea of 
mobile conduction electrons (but does not hybridize with 
them). Thus, in this example fid would represent the hole 
number operator hh ~ Wh. 

More generally, a type 2 quench presupposes a Hamil- 
tonian H{nx) that depends on a conserved charge, say fix 
[such as fi^ for the LCM], with eigenvalues [such as 
Ud — 01 1]. The Hilbert space can then be decom- 
posed into distinct, dynamically disconnected sectors, 
labelled by and governed by effective Hamiltonians 
H{nx), whose ground states have the form \G{nx)) = 
\nx)\g{nx)) ■ A type 2 quench is induced by an operator, 
say A"^ [such as (P for the LCM], whose action changes 
the conserved charge, thereby connecting two distinct 
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sectors, say {n'^\X^nx) = 1, with ^ Ux- For such 
a quench we make the identifications 

Hi = Hin,) Hi = , (9a) 

|Gi)-K)|gK)), |Gf) = K)|g«)). (9b) 

The overlap {G\\Gi) = vanishes triviaUy, because 
{nx\n'^) = 0. However, define 

|Vi)^Xt|Gi) (10) 

to be the "initial post-quench state" obtained by the ac- 
tion of the charge switching operator on the initial 
ground state. [Fig. 1(e) illustrates this state for the LCM 
with = d) ] Then the overlap 

Ox ^ m\Gi)\ = |(.9(nJ|g«))| ^ iV-^A| (n) 

again shows AO, since it is equal to the overlap of two 
Fermi sea ground states corresponding to different local 
charges. The corresponding exponent in Eq. (11) can 
again be related to a displaced charge. Ax = A^, but 
now the latter should compare the total charge within 
Marge described by the states |Gf) and IV'i): 

A|^ = (Gf|ntot|Gf) - (V'ilj^totlV'i) • (12) 

A^^ can be interpreted as the charge (in units of e) that 
fiows into Viargc during the post-quench time evolution 
from to |Gt) subsequent to the action of X'^ . To 
simplify notation, we will often omit the superscript ch 
distinguishing the displaced charge A^*^ from the AO ex- 
ponent Ax, since the two are equal in any case. 

Composite type 2 quench: Let us now consider a more 
complicated version of a type 2 quench, induced by a 
composite operator of the form = X^ . Here X^ 
switches between disconnected sectors of Hilbert space 
as above, while & does not; instead, & is assumed to 
be a local operator which acts on the dot or in the Fermi 
sea at the location of the dot, but commutes with fix- 
For the LCM, an example would be = c\ so that 
creates two electrons, one on the dot, one in the Fermi 
sea at the site of the dot. 

A physical realization hereof is furnished by the edge- 
ray edge effect occurring in X-ray absorption spec- 
troscopy (XAS), where an incident X-ray photon is ab- 
sorbed by an atom in a crystal, accompanied by the 
creation of a core hole [X^ — W) and the transfer of 
a core electron into the conduction band of the metal 
= c^)?'' Another example is the Kondo exciton dis- 
cussed in Refs. 15,16, where the absorption of a photon 
by a quantum dot is accompanied by the creation of an 
electron-hole pair on the dot, described by G^ = and 
= h\ respectively. In this example, the hole num- 
ber fih ~ h)h is conserved, but the dot electron number 
fif, — e is not, since the Hamiltonian contains dot-lead 
hybridization terms of the form (e^c-fc^e) (see Refs. 15,16 
for details). 



For a composite type 2 quench, the initial and fi- 
nal Hamiltonians and ground states are defined as in 
Eqs. (9), but the post-quench initial state is given by 

W:)^r\G,)^C^\^i). (13) 

Its overlap with the final ground state |Gf) to which it 
evolves in the long time limit has the form 

Oy EE |(V.;|G^)| = \{9{nx)\C\g'K))\ - iV"^^- . (14) 

The exponent Ay arising here is related to Ax and can 
be found using the following argument, due to Hopfield.*^ 
Due to the action of C\ the states IV'O and It/ii) describe 
different amounts of initial post-quench charge within the 
volume l^argc- We will denote the difference by 

Ac = mfitotWd - (V'ilntotl^i) . (15) 

For example, if C^ is a local electron creation or an- 
nihilation operator, then Ac = 1 or —1, respectively [as 
illustrated in Figs. 1(f) and (g)]. However, since an initial 
charge surplus or deficit at the scattering site is compen- 
sated, in the long-time limit, by charges flowing to or 
from infinity, the ground states [GJ-) and |Gf) towards 
which \ip[) and evolve, respectively, will differ only 
by one Fermi sea electron at infinity, and hence for prac- 
tical purposes describe the same local physics. In par- 
ticular, the charge within Viarge is the same for both, 
(G^|ntot|G^) = (Gf|ntot|Gf). Therefore, the total dis- 
placed charge associated with the action of Y^ is 

Ay = (G^l^totlG^) - (V'.'l^totlVA) = Ax - Ac, (16) 

where the second equality follows from Eqs. (15) and 
(12). The exponent governing the AO decay in Eq. (14) 
is thus given by Eq. (16). Since Ac is a trivally known 
integer, knowledge of Ax for a type 2 quench suffices to 
determine the AO exponents Ay for an entire family of 
related composite quenches. 

To conclude this section, we note that a type 1 quench 
can always be formulated as a type 2 quench, by intro- 
ducing an auxiliary conserved degree of freedom (say fih), 
whose only purpose is to divide the Hilbert space into two 
sectors (labelled by n/i = or 1), within which some pa- 
rameters of the Hamiltonian take two different values. 
For example, if the quench involves changing U to [/', 
this can be modelled by replacing U hy U -\- nh{U' — U) 
in the Hamiltonian. For an example, see Sec. IIIC. 

C. AO and post-quench time evolution 

After a sudden change in the local Hamiltonian, AO 
also affects the long-time limit of the subsequent time 
evolution, and hence the low-frequency behaviour of cor- 
responding spectral functions. A prominent example is 
optical absorption, for which AO leaves its im- 
print in the shape of the absorption spectrum, by re- 
ducing the probability for absorption. This is familiar 
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from the x-ray edge problem.^ In particular, in the hmit 
of absorption frequency w very close to (but above) the 
threshold for absorption, the zero-temperature absorp- 
tion spectrum has a power-law form, with an exponent 
that is influenced by AO. Recent demonstrations of this 
fact can be found in studies, both theoretical^**'^^ and 
experimental,^^ of exciton creation in quantum dots via 
optical absorption, whereby an electron is excited from a 
valence-band level to a conduction band level. 

In this subsection, we will analyse the role of AO for 
the time evolution after a type 2 quench of the form (8). 
We consider the following generic situation: For i < 0, 
a system is in the ground state \Gi) of the initial Hamil- 
tonian H[ (with ground state energy Ei), describing a 
Fermi sea under the influence of a local scattering poten- 
tial. At t = 0, a sudden change in the local potential 
occurs, described by the action the local operator . It 
switches sector to n^, yielding the post-quench initial 
state — X^Gi) at time t = 0+, and switches the 
Hamiltonian from Hi to Hf. 

The subsequent dynamics can be characterized by the 
correlator 

Gxit) = ~te'^"'9it){Gi\Xit)X^\Gi) , (17) 

where X(t) = e'^'^Xe-'^'^ reflecting the fact that X 
switches Hi to Hi. The phase factor e*""* is included 
for later convenience, with ujq to be specified below [after 
Eq.^(22)]. 

Since the Fermi sea adjusts in reaction to the sudden 
change in local potential at i = 0, AO builds up and the 
overlap function Qx{i) decreases with time. It is known 
since 1969 that in the long-time limit it decays in power- 
law fashion as^'^ 

Gx{t)^t-^i, (18) 

where Ax is the exponent governing the AO decay of 
Ox in Eq. (11). This can be understood heuristically by 
expanding Eq. (17) as 

ie-^(^'+"°)*ex(t) = 0(t)(^i|e-*^'*|Vi) (19a) 
= e{t)Y^e'^^-'mn)\\ (19c) 

n 

where \tpi{t)) = e~'^f*|?/'i) describes the time-evolution 
for t > 0, and \n) and En represent a complete set of 
eigenstates and eigenenergies of Hi. In the long-time 
limit Eq. (19c) will be dominated by the ground state 
\Gi) of Hi (with eigenenergy iJf), yielding a contribution 
KV'ilG'f)^ that scales as N'^^ [by Eq. (11)]. Now, as 
time increases, the effect of the local change in scattering 
potential is felt at increasing length scales L{t) ^ Vft, 
with Vf the Fermi velocity; regarding \Gi) as the lowest 
eigenstate of Hi in a box of size N ^ L{t), the AO of 
|(V'i|Gf)|2 _ implies Eq. (18). 



For a composite type 2 quench induced by = C^X\ 
we can conclude by analogous arguments that 

gY{t)^t~^-, (20) 

where Ay is the displaced charge of Eq. (16). 

For future reference, we also introduce the correlator 

e;gi(i) = -i9{t){G\e'"'Ce-'"'C^\G) - i^^^" (21) 

of an operator that does not switch between dynam- 
ically disconnected sectors, i.e. that commutes with n^; 
[examples of such operators are given in the discussion 
before Eq. (13) above]. Then Eq. (21) is a standard 
equilibrium correlator, with Hi — Hi, in contrast to the 
quench correlator Gx{t) of Eq. (21), where Hi ^ Hi. For 
such an equilibrium correlator the decay exponent rj"^ is 
called the scaling dimension oi . A local operator C*^ 
is relevant, marginal or irrelevant under renormalization 
if 77^ < 1, = 1 or > 1, respectively.'^^ 



D. AO and spectral functions 

Next we consider the spectral function corresponding 
to Gxit), 

Axioj) = -^3 |^Jdte'('^+*^)*fo(t)j (22a) 
= ^|(n|Xt|Gi)|25(c^-S„ + Si+c^o)- (22b) 

n 

It evidently has the form of a golden-rule transition rate 
for X^-induced transitions with excitation energy uj + ujq 
and is nonzero only for lu above the threshold frequency 
wth = (^'f — ^-i) ^ '^o- For simplicity, we will here and 
henceforth set LUth = by choosing ujq = Ei — Ei. Note 
the sum rule / dujA{uj) = {Gi\XX^Gi), which can be 
used as consistency check for numerical calculations. 

Equation (18) implies that in the limit uj — > wth = 0, 
the spectral function behaves as 

^x(^)^c.-i+2''- , 77x = iA|. (23) 

The definition of rjx is deliberately chosen such that 
Eq. (23) parallels the form of the equilibrium spectral 
function corresponding to G(j{t) of Eq. (21), namely 

A'i^iuj) ^ uj-^+^^c . (24) 

Now consider the spectral function yly(aj) involving 
the composite type 2 quench operator = C^X^. 
Equations (20) and (16) immediately lead to the pre- 
diction 

AyH - ^-1+2"- , VY = U^x ~ Ac)' , (25) 
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to be called the generalized Hopfield rule, since the 
essence of the argument by which we have obtained it 
was first formulated by Hopfield."* 

A physical situation for which Eq. (25) is relevant is the 
edge-ray edge effect occurring in X-ray absorption spec- 
troscopy (XAS). There we have = c'h^ (as explained 
above), and = 1. Thus Eq. (25) yields 

A,M ^ c^-i+(^'--i)' = ^-2A.+A,^ (26) 

reproducing a well-established result for the X-ray edge 
absorption spectrum [Ref. 4, p. 48; Ref. 5, Eq. (66)]. In 
the literature, —2Ah is often called the "Mahan contri- 
bution" to the exponent, and the AO contribution. 
Since A^, < 1, one has 2 Ah > A^, i.e. "Mahan wins", and 
Ahci^) diverges at small frequencies. For present pur- 
poses, though, it is perhaps somewhat more enlightening 
to adopt Hopficld's point of view, stated in Eq. (25), ac- 
cording to which both terms, — 2A^ and A^ arise from 
the AO exponent (A/j — l)"^. 

Equations (11), (23) and (25) will play a central role 
in this work. Their message is that the near-threshold 
behaviour of spectral functions of the type defined in 
Eq. (22) is governed by an AO exponent that can be ex- 
tracted from the overlap {ipi\Gf) between the initial post- 
quench state and the ground state \G{) to which it 
evolves in the long-time limit. 

To conclude this section, we remark that the above 
analysis generalizes straightforwardly to models involving 
several species or channels of electrons, say with index /i, 
provided that the channel index is a conserved quantum 
number (i.e. no tunnelling between channels occurs). ■^^ 
Then the initial and final ground states will be products 
of the ground states for each separate channel, so that 
Eq. (1) generalizes to 

All power laws discussed above that involve A^q (or 
quantities derived therefrom) in the exponent can be sim- 
ilarly generalized by including appropriate products over 
channels. 



E. AO exponents and NRG 

Results of the above type have been established analyt- 
ically, in the pioneering papers from 1969, Refs. 2-5, only 
for the simple yet paradigmatic case of the X-ray edge ef- 
fect. Nevertheless, Eq. (25) can be expected to hold for 
a larger class of models, as long as the setting outlined 
above applies. Indeed, it has recently been found to hold 
also in the context of the Kondo exciton.^'*"*^ The pur- 
pose of this work, therefore, is to establish the validity 
of the connections between Eqs. (11), (23) and (25) for a 
series of models of increasing complexity. We shall do so 
numerically using NRG, since for most of these models 



an analytical treatment along the lines of Refs. 2 and 5 
would be exceedingly tedious, if not impossible. How- 
ever, the requisite numerical tools are available within 
NRG,^^''^° and have become very accurate quantitatively 
due to recent methodological refinements. ^^''^^''^^ 

NRG, developed in the context of quantum impurity 
models, offers a very direct way of evaluating the overlap, 
since it allows both ground states \Gi) and \Gf) to be 
calculated explicitly. Models treatable by NRG have the 
generic form H — Hb + Hd- Here 

describes a free Fermi sea involving channels of 
fermions, with constant density of states p per channel 
and half-bandwidth D = l/(2p). (When representing 
numerical results, energies will be measured in units of 
half-bandwidth by setting D = 1.) Hd, which may in- 
volve interactions, describes local degrees of freedom and 
their coupling to the Fermi sea. 

Wilson discretized the spectrum of Hq on a logarithmic 
grid of energies ±DK-^ (with A > 1, /c = 0,1,2,...), 
thereby obtaining exponentially high resolution of low- 
energy excitations. He then mapped the impurity model 
onto a semi-infinite "Wilson tight-binding chain" of sites 
fc = to oo, with the impurity degrees of freedom coupled 
only to site 0. To this end, he made a basis transforma- 
tion from the set of Fermi sea operators {c^^} to a new 

set {/fc^}, with /o^ oc = ?/'m(0) = Y.e ^em' chosen such 
that they bring Hq into the tridiagonal form 

Tic OO 

^B^EE*'^-(4A-i,M+^-'=-), (29) 

^l=l k=i 

with hopping matrix elements tk oc DK~^/^ that decrease 
exponentially with site index k along the chain. Because 
of this separation of energy scales, the Hamiltonian can 
be diagonalized iteratively by solving a Wilson chain of 
length k (restricting the sum in Eq. (29) to the first k 
terms) and increasing k one site at a time. The number 
of kept states at each iteration will be denoted by . 

For a Wilson chain of length k, the effective level spac- 
ing of its lowest-lying energy levels is set by the small- 
est hopping matrix element of the chain, namely A"*^/^; 
such a Wilson chain thus represents a real space system 
of volume Marge ^ A*^/^. Thus, the overlap between the 
two ground states of a Wilson chain of length k can be 
expressed as^^ 

|fc(Gi|Gf)fe|~A-t'^Ao =e-"^ (30) 

where a = (logA/4)A^Q. Explicit calculations show^^ 
that an exponential decay of the form Eq. (30) applies for 
the overlap between any two states \E\)-k and li^f)^; rep- 
resenting low-lying excitations w.r.t. |Gi)fc and |Gf)fc at 
iteration k, respectively. More technically, k(J^\\E'i)k ~ 
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g-afc jj^Qjfjg whenever \Ei)k and \E{)k represent NRG 
eigenstates with matching quantum numbers from the 
k-th NRG shell for H[ and Hf, respectively, and their 
overlap is calculated for increasing k. For multi-chain 
models, we note that channel-specific exponents such as 
^AO,/i [see Eq. (27)] can be calculated, if needed, by con- 
sidering Wilson chains with channel-dependent lengths. ■^^ 

Within the framework of NRG, a consistency check is 
available for the value of Aao extracted from Eq. (30): 
Aao should be equal to the displaced charge Ach of 
Eq. (7), which can also be calculated directly from NRG 
by calculating the expectation value of ntot for |Gi) 
and \Gf) individually.^^ This check was successfully per- 
formed, for example, in Refs. 14 and 15, within the con- 
text of the single impurity Anderson model; for a recent 
systematic study, see Ref. 21. We have also performed 
this check in the present work wherever it was feasible. 

Within NRG, it is also possible to directly calculate 
spectral functions such as Axi^) of Eq. (22). To this 
end, one uses two separate NRG runs to calculate the 
ground state |Gi) of Hi and an approximate but com- 
plete set of eigenstates |n) of H{.^^'^^ The Lehmann sum 
in Eq. (22) can then be evaluated explicitly, while 
representing the (5-functions occurring therein using a log- 
Gaussian broadening scheme. To this end, we follow the 
approach of Ref. 34, which involves a broadening param- 
eter a. (The specific choice of NRG parameters A, Nk 
and a used for spectral data shown below will be spec- 
ified in the legends of the corresponding figures.) That 
this approach is capable of yielding spectral functions 
whose asymptotic behaviour shows power-law behaviour 
characteristic of AO has been demonstrated recently in 
the context of the Kondo exciton problem. ^''"^^ In the ex- 
amples to be discussed below, we will compare the power- 
law exponents extracted from the asymptotic behaviour 
of such spectral functions to the values expected from 
AO, thus checking relations such as Eq. (23) for Ax{^) 
and Eq. (25) for yl.y(w). 



III. INTERACTING RESONANT LEVEL 
MODEL 
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Figure 2: (Color online) Numerical results for the LCM of 
Eq. (2), for the type 2 quench of Eq. (8), whose initial, final 
and post-quench initial states \Gi}, |Gf) and are depicted 
in Figs. l(b,c), 1(d) and l(e-g), respectively, (a) Comparison 
of the decay exponent Ad obtained from Eq. (31) (crosses) 
with the displaced charge A^*" from Eq. (32) (pluses), for a 
number of different values of U. The two values agree very 
well (they differ by less than 0.1%), also with the analytic 
prediction Eq. (33) (solid line). As expected, A^j — >■ —1/2 
for U — >■ 00. (b) Comparison of two ways of determining 
the AO exponents rj that govern the low-energy asymptotic 
behaviour A ~ oj"^"'"^'' of the spectral functions of Eqs. (35), 
related to Figs, l(e-g): exponents obtained by fitting a power 
law to the corresponding spectra [shown in (c)] are shown 
as crosses (marked "spec", for "spectra"); the corresponding 
exponents expected from Eq. (35), using the results of (a) for 
Ad, are shown as dots (marked "exp" for "expected"). We 
find a maximal deviation of less than 1%. Here and in all 
similar figures below, the dashed lines are only guides to the 
eye. (c) Asymptotic low-frequency dependence of the spectra 
Eqs. (35), for (7 = 1, on a double logarithmic plot, allowing 
the corresponding exponents rj to be extracted. 



In this section we consider the effect of AO on dy- 
namical quantities in the context of the spinless inter- 
acting resonant level model (IRLM).^^''^^ (The effects of 
AO for some static properties of this model were studied 
in Ref. 36.) The purpose of this exercise is to illustrate 
several effects that will be found to arise also for more 
complex models considered in subsequent sections. The 
IRLM involves a single localized level, to be called dot 
level (alluding to localized levels in a quantum dot), in- 
teracting with and tunnel-coupled to a single Fermi sea. 
We consider first the case without tunnelling, in which 
case the IRLM reduces to the LCM introduced in Sec. II 
above, where adding an electron to the dot at time t = 
constitutes a type 2 quench. This leads to AO between 
the initial and final ground states, and corresponding 



nontrivial AO power laws, oj"^"'"^'', in spectral functions. 
We then turn on tunnelling, which connects the sectors 
of Hilbert space for which the dot is empty or filled, and 
hence counteracts AO. Correspondingly, the power-laws 
get modified at frequencies smaller than the renormal- 
ized level width, lo < Fion, where the AO behaviour is 
replaced by simple Fermi liquid behaviour; the effects of 
AO do survive, however, in a regime of intermediate fre- 
quencies, Fj-on < < D. Finally, we consider quenches of 
the position of the dot level, in which case AO reemerges. 
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A. Without tunnelling: LCM 

In this subsection we present numerical results for 
the IRLM without tunnelling, corresponding to the lo- 
cal charge model of Eq. (2), depicted in Fig. 1(a). We 
consider the type 2 quench of Eq. (8), with = d) . 
The initial and final ground states \Gi) and |Gt) are il- 
lustrated in Figs. l(b,c) and 1(d), respectively, and the 
post-quench initial state j-i/'i) = d^Gii in Fig. 1(e). With 
these choices the overlap K^/iijCf)! of Eq. (11) becomes 

^ |(GoMlGi)| = |(5o|5i)l ^ N-'^"^'^ . (31) 

The corresponding displaced charge obtained from 
Eq. (12) is 

= (.9i|n-sca|gi) - (.go|"-sca|5o) , (32) 

since \G{) and 1-01) describe the same dot charge, — 1. 

We used NRG to calculate the overlap Od of Eq. (31) 
and extract the exponent A^; from its exponential decay 
with Wilson chain length [Eq. (30)], for several values of 
U. As consistency check, we also calculated the displaced 
charge /S.'f- [Eq. (32)]. As shown in Fig. 2(a), the results 
for Ad (crosses) and A^*^ (pluses) agree very well. The 
displaced charge AJ^*^ is < 0, since the repulsive interac- 
tion pushes charge away from the local site. Its magni- 
tude [A^^l depends on the interaction strength: as U is 
increased from to oo, the displaced charge goes from 
to — i , reflecting the complete depletion of the initially 
half-filled Wilson chain site directly adjacent to the dot 
site [compare Figs. 1(b) and 1(d)]. Figure 2(a) shows 
that the numerical results for A^ and AJ^*^ (symbols) also 
agree with the analytical result (solid line) obtained for 
the phase shift obtained from elementary scattering the- 
ory [see e.g. Ref. 11, Eq. (25.29)], 

Arf = --tan~i {-KpU) , (33) 

with p the density of states in the Fermi sea (cf. Sec. II E). 

To study the influence of AO on dynamical quantities, 
we consider simple and composite type 2 quenches in- 
duced by acting on the initial ground state |Gi) — [Go) 
with the operators 

= , = c^d) , = cd) . (34) 

All three operators describe transitions between the — 
and 1 sectors. The analysis of Sec. II D applies directly, 
with the identifications Hi = Hq and H{ = Hi, while 
Ac = ±1 for Y^ or Y2 ^ respectively [see Figs, l(e-g)]. In 
particular, Eqs. (23) and (25) imply: 

Ad{uj) ^ ^^-1+2"^ 7?, = iA2 , (35a) 

^de(^) ^ r,d, = i(Ad-l)2, (35b) 

Adc^iuj) ^ c^-i+^^-t r,d,, = i(Ad + I f . (35c) 

Using NRG, we calculated these three spectra for several 
values of U (cf. Fig. 2(c)). In the limit of small uj, the 



(a) 




Figure 3: (Color online) (a) Cartoon of the Hamiltonian 
(36) for the IRLM. (b) The renormalized level width Frcn, 
calculated via the dot's charge susceptibility, I/tyxc^^ (dots) 
or via Eq. (37) (solid line), shown as a function of U for = 
and several values of F. As U increases from 0, Fren/F begins 
to differ significantly from its initial value, namely 1, only 
once U becomes comparable to the band-width, reaching its 
maximal value (F/D)"^/^ for U D. 

spectra show clear power law behaviour, The 
exponents rj^, rj^c, ildct extracted from these spectra are 
shown in Fig. 2(b) (crosses, marked "spec", for "spec- 
tra"). They agree well with the values expected (dots, 
marked "exp", for "expected") from Eqs. (35), based on 
the value for extracted from Eq. (31). Thus, all ways 
of determining are completely consistent, confirming 
the validity of the above analysis. 



B. With tunnelling: IRLM 

The previous subsection focused on a switch between 
two sectors of the Hilbert space, with rid — and Ud — 1, 
that were not coupled dynamically, but governed instead 
by two distinct Hamiltonians, Hi and H{. In the present 
subsection, we consider the case that the sectors with 
Ud — and Ud = 1 are coupled by tunnelling between dot 
and lead, so that the notion of an initial and final Hamil- 
tonian, acting in decoupled sectors of Hilbert space, does 
not apply. The dynamics is governed instead by the single 
Hamiltonian Hi = Hf = -ffiRLM, given by [see Fig. 3(a)] 

^IRLM = SdM +U{M - 1/2) (c^c -1/2) (36) 

+ ^ ""^eCe + + c\d) . 

We assume, here and in all later settings, that the hy- 
bridization of the dot level with the Fermi sea states is e- 
independent, with F being the bare width of the dot level. 
Here, in contrast to the local charge model of Eq. (2), the 
interaction term is taken to be particle-hole symmetric, 
so that the model is particle-hole symmetric for Ed = 0. 
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The presence of the interaction, U, is known to effec- 
tively modify the level width, both by reducing the 
density of states of the leads near the dot, and by induc- 
ing AO in the leads when the dot occupancy changes. 
The precise interplay between these effects can be quite 
intricate and was studied in Ref. 36. A quantitative anal- 
ysis can be performed by defining a renormalized level 
width in terms of the charge susceptibility, Ficn = I/ttXc- 
At the point of particle-hole (ph) symmetry (e^ = 0), an 
analytic formula for the latter is available, '^^ 



r,,jD = (T/D) 
where is given by 



l/(2-(l+Af 



Af = --tan- 

TT 



(37) 



(38) 



can be interpreted as the change in scattering phase 
shift that a system with F = 0, Ed — experiences if the 
local occupancy is changed abruptly from n^j = to 1. 
The form of Eq. (38) is analogous to Eq. (33) for A^, 
with two differences: since the final scattering potentials 
being compared have amplitude —U/2 and U/2 (instead 
of and U), the argument of tan~^ has an extra factor 
of 1/2, and there is an extra prefactor of 2. 

The dependence of Fren on U is illustrated in Fig. 3(b), 
which shows good agreement between the NRG results 
for 1/nXc (dots) and the analytic formula (37) (fines). 
For U much smaller than the bandwidth D, Fj-en/F is es- 
sentially equal to 1; it strongly increases once U becomes 
of the order D, and saturates to {T/Dy-/'^ for U :$> D. 

Let us now consider the equilibrium spectral functions 
for the operators of Eq. (34), A^"^, A2 and They 
are defined as in Eq. (22) but with Hf = Hi, because for 
the IRLM, where Ud is not conserved, none of these oper- 
ators induces a quench. Therefore, the cj — ^ behaviour 
of their correlators is expected (and indeed found) to be 
independent of AO. However, quite remarkably, traces 
of AO do show up in an intermediate frequency regime, 
uj* < u! < D, where 1/uj* corresponds to the time scale 
within which charge equilibration takes place. Below the 
energy scale w* the quantum impurity becomes strongly 
correlated with the Fermi sea and for the present model 
we have uj* ~ Frcn- Let us therefore discuss the two 
regimes, uj below or above a;*, separately. 

In the regime a; ^ w*, the spectral functions are found 
to have the following asymptotic form A'^'^ ~ ui' 
[cf. Figs. 4(a-c)]: 







cq 

Vd 


= 1/2, 


(39a) 






cq 
Vdc 


= 2, 


(39b) 






cq 


= 2 . 


(39c) 



The exponents arising here can be understood analyti- 
cally using elementary, though not entirely trivial argu- 
ments, based on the fact that the lowest-lying excitations 
of this model have Fermi liquid properties. We refer the 
reader to the Appendix for a detailed analysis. 
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Figure 4: (Color online) (a)-(c) The equilibrium spectral 
functions A^iuj), AHii^) and A2i {^) for the IRLM, showing 
a crossover from trivial power laws, uj^^'^^'^ '', for cj < oj*, to 
AO power laws, uj~^'^^^ , for lj* < lj < D, with the crossover 
frequency uj* given by Frcn- (d) Comparison of the exponents 
ri"'^ (triangles) and 77" or rj'^ (crosses) extracted from the spec- 
tra shown in (a-c), with the values expected from Eqs. (39) for 
rj°^ (solid lines), and from Eqs. (41) for 77° or from Eq. (40) 
for rj'j (dashed lines with dots), for several values of U. In 
(a), arrows indicate the scale cj* that separates the regimes 
UJ* < UJ < UJ* and uj* < u < D, where A^ scales according to 
Eqs. (40) or (41a), respectively. 



Now consider the regime lo* < uj < D. As shown in the 
corresponding regime of w/Fj-cn > 1 h^ Figs. 4(a-c), each 
of the equilibrium spectral functions A^, A2 and A^f^ 
exhibits another, different power-law there. For A^ we 
actually find that within this regime two different power- 
laws can be discerned: First, in a regime uj* < uj < Q* 
we find, 
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where A^*^ is given by Eq. (38). The exponent i]'^ cor- 
responds to the leading correction for weak interactions 
{U/D <C 1) to the pure Lorentzian decay of the spectral 
function of the d-level (see Eq. (A. 2) of the Appendix), 
as can be shown using methods discussed in Refs. 37 and 
38. The scale ui* that sets the upper limit for this be- 
haviour is marked by arrows in Fig. 4(a) and decreases 
with increasing U /D. For U/D sufficiently small that w* 
lies far below the bandwidth D, we find a second power 
law within the window Q* < uj < D, namely 

For the other two spectral functions we find throughout 
the regime lu* < uj < D: 

ATM - Vl = ^(Af - ir , (41b) 

^^,^,(^)^..-i+<t, r^l.^^-iAf + lf. (41c) 

Remarkably, Eqs. (41) have the same form as Eqs. (35), 
except that is replaced by A^'' of Eq. (38), i.e. by 
the AO exponent involved in abruptly changing the lo- 
cal occupancy from to 1 (in the absence of tunnelling) . 
That this exponent should emerge is natural, since the 
corresponding correlators Gd, Qdc and all involve an 
operator d) that places an electron on the dot at time 
t = Q. Although the dot occupancy nd{t) will relax 
back to its initial value nj; in the long time limit, this 
requires times t ^ 1/w*. In contrast, the lead electrons 
react to the change in local charge on the much shorter 
time scale \/D. Thus, in the window of intermediate 
times, 1/D <C t ^ corresponding to frequencies 

w* <C w ^ I?, the situation is similar to that of the pre- 
vious subsection, where we had F = and a change in 
dot occupancy from to 1 induced changes in the lead 
phase shifts, accompanied by AO. Thus, the exponents 
rp arising in Eq. (41) can be identified as the (equilib- 
rium) scaling dimensions of the corresponding operators 
calculated in the absence of tunnelling (which is why we 
use a superscript on such exponents, here and below). 
This explains the similarity between the behaviour de- 
scribed by Eqs. (41) and Eqs. (35). Note that the scaling 
dimension 77^^^ [Eq. (41c)] of the tunnelling operators dc^ 
and &d satisfy < rf^^^ < 1/2 [since for [/ > 0, we have 

— 1 < IsP^ < 0, by Eq. (38)], thus tunnelling is always 
relevant for this model. 

We conclude this subsection with a comment on the 
fact that A'T{uj) crosses over from non-AO behaviour 
[Eq. (40)] to AO behaviour [Eq. (41a)] as U/D is in- 
creased past 1. AO behaviour is absent for U/D <^ 1 
because this situation corresponds essentially to a non- 
interacting resonant-level model, for which Gd'^{t) does 
not show power-law behaviour of the type assumed in 
Eq. (18); instead it decays exponentially (^ e^^*), caus- 
ing the spectral function A'T{lu) to have an essentially 
Lorentzian form. Equation (40) is the large-frequency 
limit of the latter, but including the leading corrections 
in U/D, calculated using methods discussed in Refs. 37 
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Figure 5: (Color online) (a) Cartoon of the quench which 
occurs when an electron- hole pair is created at time t = 0, 
see Eq. (43). (The cartoon depicts the situation relevant for 
exciton creation by absorption of a photon, which excited an 
electron from a valence-band to a conduction band level of a 
semiconducting quantum dot.) (b) The exponent Ah [from 
Eq. (11)] and the displaced charge A^*" [from Eq. (12)], for 
the quench of Eq. (42), as function of the quench range W. 
(c) Corresponding values of the AO exponents r/hd, Vhdc and 
J7h^j,t, extracted from the asymptotic behaviour u"^"'"^'' of 
spectral functions (crosses), or as expected from Eqs. (46) 
(dots). Typically, relative errors are less than 1%. 



and 38. However, once U/D becomes > 1, AO does be- 
gin to matter, implying a regime of power-law decay for 
Gd^{t) on intermediate time scales, leading to Eq. (41a) 
for ^7(0.). 



C. Quantum quench of level position 

In the previous subsection we emphasized the impor- 
tance of the scale uj* , which separates the low- and 
intermediate- frequency regimes, showing trivial expo- 
nents or AO exponents, respectively. It is instructive 
to study the role of the scale ui* in a slightly different 
but related context, namely quench spectral functions 
involving a quantum quench of the level position. This 
will shed further light on the AO between states with 
different local level occupancies. 

Concretely, we consider initial and final Hamiltonians 
that both are of the form Eq. (36), but with initial and 
final level positions that are symmetrically spaced on op- 
posite sides of the model's symmetry point at Sd = 0: 

e'd^W/2 '^^'^ e^^ = -W/2. (42) 
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Figure 6: (Color online) (a) The AO exponent Ah [extracted 
according to Eq. (11)] as function of VK/Frcn, for several 
values of U. For IF/Fron 3> 1, Ah approaches its maxi- 
mal value A™^^. As W is reduced below Fren, A^ drops 
below its maximal value and decreases with W [linearly so 
for W/Tren <C 1] • (b) The final and initial occupancies 
(squares, upper curves) and n\ (circles, lower curves) as func- 
tions of W/T-ra-n, for the same values of U [same color code 
as in (a)], (c) The maximal value A™^ of the AO exponent 
Ah, extracted from the W/Fren 1 regime of (a) (crosses, 
"calc"), or expected from Eq. (44) (solid line, "exp"); the 
relative deviations are well below 1%. 



Although this is an example of a type 1 quench, it will be 
convenient (mainly for notational reasons) to reformulate 
this situation as a type 2 quench. To this end we use the 
Hamiltonian 

H = W {I /2 - nh)nd + U [fid- I /2){c^c- 1/2) 




(43) 



where we have introduced an auxiliary degree of freedom, 
called "hole" (in analogy to the role of holes in exciton 
creation by optical absorption^^-^®), with hole counting 
operator fih = Wh. The hole has no dynamics; its only 
role is to distinguish two distinct sectors of Hilbert space, 
in which the dynamics is described by Hi or iff, with hole 
number rih — or 1, respectively [see Fig. 5(a)]. The type 
2 quench that switches between these sectors is induced 
by = . The overlap Oh ~ N~2^h between the 
initial and final ground states is characterized by an AO 
exponent A^, [Eq. (11)] that is equal to the charge Af^ 
displaced by the quench [Eq. (12)]. 

The magnitude of A/j increases with the range W of the 
quench, as shown in Fig. 5(b) (linear scale) and Fig. 6(a) 
(log-log scale). Note, in particular, that the scale on 



which the quenching range, needs to change in order 
for the AO exponents to change significantly, is given by 
Fron- This is natural: when W ^ Frcn, the two states 
\Gi) and \Gf) connected by the quench describe dots with 
strongly different occupancies, njj ~ vs. ~ 1, see 
Fig. 6(b). Hence the AO [Eq. (11)] of the corresponding 
Fermi seas will be strong. The maximum possible value 
of the exponent Ah is 



(44) 



with Af'iU) < given by Eq. (38). The first term simply 
gives the J7 — >■ oo value of the change in dot occupancy 
induced by the quench, namely 1; the second term reflects 
the reaction of the Fermi sea to this change, cf. Sec. IIIB. 

Following the arguments of Sec. II D, the nonequilib- 
rium spectral functions ^y(a;), defined for 



= cd^^h) , (45) 



are expected to show the following AO behaviour for 
a; ^ 0: 

AdH ^ ^.-1+2"- r^hd^\{Ah-l)\ (46a) 
AdcM ^ 0."'+""'''^ flhdc=\{Ah-2)\ (46b) 
Adct {lo) ^ c.-i+^^.^ct ^^^^t - \Al . (46c) 

The reason for the specific form of the exponents is that 
for the correlators Qhdi Ghdc or GhdcU at t = the local 
charge (on the d-level or in the Fermi sea) is increased 
by one, two or zero, respectively [i.e. Ac = 1, 2 or in 
Eq. (16)]. Figure 5(c) shows that the exponents (crosses) 
extracted from the asymptotic behaviour Ayd^) are in- 
deed in good agreement with values expected (dots) from 
Eqs. (46). 



IV. POPULATION SWITCHING WITHOUT 
SENSOR 

The models investigated so far served as testing ground 
for the influence of AO on various types of spectral func- 
tions. The following two sections have the concrete mo- 
tivation to clarify the role of AO in the context of quan- 
tum dot models that display the phenomenon of popula- 
tion switching (ps). 17,18, 39-44 such models, a quantum 
dot, tunnel-coupled to leads, contains levels of different 
widths and is capacitively coupled to a gate voltage that 
shifts the levels energy relative to the Fermi level of the 
leads. Under suitable conditions, an (adiabatic) sweep of 
the gate voltage induces an inversion in the population 
of these levels (a so-called population switch), implying 
a change in the local potential seen by the Fermi seas in 
the leads. Goldstein, Berkovits and Gefen (GBG) have 
argued in Ref. 17,18 that in this context AO can play an 
important role. In particular, they pointed out that for 
a model involving a third lead acting as a charge sensor, 
the effects of AO can be enhanced to such an extent that 
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population switching becomes abrupt, i.e. turns into a 
phase transition. Our goal is to elucidate the influence 
of AO by using the tools developed above in the context 
of the IRLM. 

In the present section, we will study population switch- 
ing in a two-lead model (without charge sensor) , which is 
equivalent to an anisotropic Kondo model. The 
corresponding Kondo temperature, Tk, sets the width of 
the population switch as function of gate voltage. We 
calculate the spectral function A'p{u!) of the pseudospin- 
flip operator and show that Tk also acts as the crossover 
scale Lo* that separates a low-frequency regime showing 
Fermi-liquid power laws from an intermediate-frequency 
regime revealing AO exponents. We investigate the ori- 
gin of the latter by a quantum quench analysis similar 
to that of Sec. IIIC above. In the following section, we 
will generalize the model by adding a charge sensor and 
analyse how this enhances the effects of AO. 



A. Width of switching regime 

We consider a model involving two single-level dots 
(/i = L, R) and for convenience choose their level ener- 
gies Ed to be equal, so that the PS always occurs at the 
particle-hole symmetric point, = —U/2. (Note that PS 
occurs also for nondegenerate levels, as long as their level 
spacing is smaller than the difference of their level widths 
I^.) The levels have an electrostatic coupling U > and 
are each tunnel-coupled to its own lead [see Fig. 7(a)]: 



(a) 



(47) 



SfJ. /J V £ 

(We use notation analogous to that of Sec. IIIB.) We 
choose the level widths to be strongly asymmetric and 
will use a fixed value of their ratio, Fr/Fl = 20, through- 
out. The model thus has the form of a spin-asymmetric 
single- impurity Anderson model (SIAM), where /i acts as 
pseudospin index. 

As illustrated in Fig. 7(b), this model shows PS when 
Ed is decreased past Ed = —U/2 (the particle-hole sym- 
metric point): as this "switching point" is crossed, the 
occupancy of the broad level (dashed lines) changes from 
near 1 to near 0, and vice versa for the narrow level (solid 
lines). We define the width of the switching regime, Wps, 
as the difference, 
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Figure 7: (Color online) (a) Cartoon of the Hamiltonian (47) 
for the asymmetric SIAM. (b) The occupations til (solid lines) 
and nR (dashed lines) of the left and right level, respectively, 
as functions of £d, for several values of Fa, at a fixed ratio of 
Fr/Fl = 20. As Ed is lowered past the particle-hole symmet- 
ric point at Ed ~ —U/2, population switching occurs, with ur 
changing from near 1 to near 0, and vice versa for til. Inset: 
zoom into the switching region around Ed = —U/2, showing 
that the switch is continuous (as function of Ed) even though 
the switching region becomes narrower for decreasing Fr. (c) 
Comparison of fcoIFps [from Eqs. (48) and (51), crosses], Tk 
[from Eq. (52), solid], and the inverse pseudospin suscepti- 
bility 1/Xs (pluses). All three quantities evidently decrease 
similarly with decreasing Tr/U. 



Figure 7(b) and its inset show that the width of the 
switching regime decreases with decreasing T^, without, 
however, dropping to zero as long as the level widths 
are nonzero. This behaviour can be understood as 
follows. In the vicinity of the particle-hole sym- 
metric point, only two local charge configurations are rel- 
evant, namely those with occupancies (wli'^r) equal to 
(0, 1) or (1,0). The spin-asymmetric SIAM can thus be 
mapped onto an anisotropic Kondo model by a Schrieffer- 
Wolff transformation. This leads to an anisotropic pseu- 
dospin exchange interaction of the form 



between those two values of Ed, located symmetrically on 
either side of the switching point, at which the occupation 
of the right level is nR+ = ^n'^^ (> i) or nR_ = 1 — 
nR+ (< 5), respectively, where n^'^ is the largest value 
reached by riR for Ed > —U/2, to the right of the PS. 



-ffcxch = (?^l + ?/r)(44-^^r^r)(c[cl-4cr) 
+2VCW^(cl444 + li-c-) 
+i?cff(4'^"L-44)/2, (49) 
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respectively, with coupling constants given by 
T„ f I l_ 



TT 



and effective magnetic field 

Soff = 6o(ed + C//2), ha 



4(r, 



R - -^Lj 



(50) 



(51) 



The corresponding Kondo temperature is given by the 
following expression:^^ 



^ v/t/(rL + Tr) 

Tk = exp 



7r£rf(£/ + £rf) ^^Tl 

2i7(rL-rR) ""tr 



Note that Tk decreases exponentially if IJ; is decreased 
with a fixed ratio of Fr/Fl and actually becomes zero 
for 1^ = (the argument of the exponent in Eq. (52) is 
negative, since < 0). 

Now, Tk can be associated with the energy gained 
by forming a ground state involving a screened local 
pseudospin, which in the present setting translates to a 
ground state involving a coherent superposition of config- 
urations with local occupancies (0, 1) and (1,0). Screen- 
ing will cease when Sd deviates sufficiently from the sym- 
metry point —U/2 that the effective magnetic field |i?cff| 
exceeds Tk, in which case the ground state will be dom- 
inated solely by the (0, 1) or (1, 0) configuration, instead 
of involving a coherent superposition of both. Thus the 
switching width will be set by &oW^ps — up to a nu- 
merical constant of order unity. 

Figure 7(c) confirms this expectation. It shows that 6o 
times the switching width Wps [from Eq. (48)] (crosses) 
and the Kondo temperature Tk at e;^ = — J7/2 [from 
Eq. (52)] (sohd line), when plotted as functions of Fr/C/ 
at fixed Fr/Fl, are indeed almost perfectly proportional 
to each other. As a numerical consistency check, Fig. 7(c) 
also shows the inverse of the zero-temperature pseu- 
dospin susceptibility of the dot levels, 1/xs (pluses), con- 
firming that Tk = 1/Xs- (This is analogous to the rela- 
tion Fj-on = ^/t^Xc of Sec. IIIB.) 



B. AO in dynamics of pseudospin-flip operator 

Let us now explore the role of AO in population switch- 
ing. To this end, we note that the effective exchange 
interaction i^oxch of Eq. (49) is similar in structure to 
the IRLM of Eq. (36) : both involve two charge configu- 
rations (the former (0,1) and (1,0), the latter and 1), 
which induce different phase shifts in the leads due to a 
dot-lead interaction term, and which are connected by a 
tunnelling term. More formally, the relation between the 
IRLM and PS is revealed by the equivalence of both mod- 
els to the Kondo model (for the IRLM, this equivalence 
is discussed, e.g., in Refs. 11,35,36). Thus, we may ex- 
pect AO to play a similar role for both models, and hence 
perform an analysis similar to that in Sections III B and 

inc. 
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Figure 8: (Color online) (a) The pseudospin-flip spectral 
function (cj) [cf. Eq. (53)] for the PS model without charge 
sensor, for several values of T-n/U with fixed ratio Fr/Fl, 
calculated at Sd — —U/2: plotting A'^{uj)/A'^{Tk) versus 
uj/Tk yields a scaling collapse. The frequency dependence of 
the curves qualitatively changes at Tk'- for lj < Tk we find 
Fermi liquid behaviour, ~ uj"^ (dotted line), while for lj > Tk 
each curve shows a nontrivial AO power-law, ~ uj~^'^^^y 
[Eq. (54)], exemplified by the dashed and dash-dotted lines 
for Fa/iy = 0.2 and 0.05, respectively, (b) Comparison of 
the values for riy expected from Eq. (54b) (solid line), or ex- 
tracted from the spectral function ^^(tj) by using Eq. (54a) 
in the intermediate-frequency regime between Tk and the high 
frequency maximum (crosses). The accuracy of this extrac- 
tion deteriorates with increasing Fr/?/, since this reduces the 
width of the latter regime [see (a)], causing the relative er- 
ror between crosses and solid line to increases from 1% for 
Fr/17 = 0.05 to 10% for Fn/U = 0.2. 



Specifically, let us study the spectral function (w) 
of the pseudospin-flip operators occurring in iJcxch, 



(53) 



These induce transitions between the local charge con- 
figurations (0,1) and (1,0) and simultaneously add an 
electron to one lead while removing an electron from 
the other. (Such a transition does not constitute a 
quench, since for the present model is not conserved.) 
A'p should, in some respects, be analogous to A^cf of 
Sec. IIIB. We have thus calculated A'P numerically, us- 
ing the Hamiltonian -ffsiAM of Eq. (47). Indeed, Fig. 8(a), 
which shows A'^{uj) for several values of Fr/C/, exhibits 
several features reminiscent of Fig. 4(c) for A^cii^)'- (i) 
A crossover scale cj*, separating a regime of very low fre- 
quencies from one of intermediate frequencies, is clearly 
discernible; it is given by lj* ~ Tk. (ii) When prop- 
erly rescaled by plotting A'P{uj)/A^{u}*) versus uj/uj*, 
all curves collapse onto each other in the regime oj < D. 
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(iii) In the low-frequency regime u uj* we find the 
same Fermi-Uquid power law, A'^{ui) ^ (dotted line), 
as for y^^^i (w) [cf. Eq. (39c)]. (An analytical explana- 
tion for this fact is given in at the end of the Appendix.) 

(iv) In an intermediate- frequency regime i^* ^ w ^ Whc, 
whose upper limit cjhc is a high-energy scale set by the 
minimum of the bandwidth or the cost of local charge 
fluctuations, we find an AO-dominated power law, 

A^^{uj) - cj-i+2<, cj*<w<cjhc. (54a) 

Though the numerical calculation of A'^{lo) was per- 
formed using the full Hamiltonian /fsiAM of Eq. (47), 
tunnelling is not important on the short time-scales that 
govern the frequency regime lo > uj* . Hence, we expect 
the exponent rf^ found from Eq. (54a) to be equal in 
value to that which one would obtain in the w — > limit 
of a calculation performed in the absence of pseudospin- 
flips, i.e. using i?exch without the pseudospin-flip terms 
in the second line of Eq. (49) (but retaining its first and 
third lines). 

Without pseudospin-flips, the correlator involving Y 
would actually constitute a type 2 quench correlator, 
because changes (riL — ?^R), which is a conserved 
quantum number for Hamiltonians without pseudospin- 
flips. Therefore, the expected value of rjy can be pre- 
dicted using the generahzed Hopfield rule [Eq. (25)]. For 
the present case of two channels that are not intercon- 
nected by tunnelling, so that the total charge within each 
channel is conserved, it can be applied to each chan- 
nel separately, adding the corresponding exponents'^ ^ [cf. 
Eq. (27)]: 

< = i(AL + l)2 + i(-AR-l)2. (54b) 

Here describes the change in phase shift, divided by 
TT, induced in lead /i by a pseudospin-flip; it is given by 
Eq. (38), with U replaced by [/^ [from Eq. (50)]. The ap- 
plicability of these arguments is confirmed by Fig. 8(b), 
which shows that the exponents extracted from the nu- 
merical spectra (crosses) agree quite well with the values 
expected from Eq. (54b) (solid line). The quality of the 
agreement deteriorates with increasing Fr/C/, because 
this reduces the width of the intermediate-frequency 
regime, making an accurate extraction of rf^ horn A'^{uj) 
increasingly difficult. 

Equation (54b) allows us to understand why PS is al- 
ways continuous in this model: Since — 1 < A^ < 0, 
the scaling dimension of Y satisfies ?7y < 1, implying 
that this operator always remains a relevant perturba- 
tion that does not flow to zero at low energy scales. This 
means that AO, although present, is not strong enough 
to completely suppress the amplitude for pseudospin-flip 
transitions. Hence, the two sectors (0,1) and (1,0) are 
always coupled by the effective low-energy Hamiltonian, 
so that population switching is continuous. 




W/ 



Figure 9: (Color online) AO for the PS model without charge 
sensing. The exponent Aao [extracted from Eq. (1)] (solid 
lines with dots) is shown as function of quench size W/U 
[Eq. (55)], for several values of Tyi/U with fixed ratio Fr/Fl, 
showing that AO becomes significant once W increases past 
TK/bo (indicated by dashed vertical lines). The exponent 
Aao increases linearly with W for W ^ T^/bo, and satu- 
rates to a maximal value of \/2 [Eq. (56)] for W ^ Tk/So. 
The corresponding values of Ach [from Eq. (7)] are not shown 
but agree with Aao with relative errors of a few percent. 

C. AO induced by quench of level positions 

As mentioned above, the operators Y^ and Y connect 
two configurations with different local occupancies, (0,1) 
and (1,0). To shed further light on the AO between such 
configurations, we now perform a quantum quench anal- 
ysis similar to that of Sec. IIIC. We consider a type 1 
quench, iJj — )• Hf, induced by changing the level posi- 
tion Ed from a value above the symmetry point, favouring 
(0,1), to one below, favouring (1,0): 

e^ = -U/2 + W/2 '^^'^ el = ~U/2-W/2. (55) 

The corresponding ground states, \Gi) and [Gf), will dis- 
play AO as in Eq. (1). Based on the lessons learnt from 
Sec. IIIC, the corresponding exponent Aao will increase 
with the width W of the quench. Indeed, Fig. 9 [to be 
compared with Fig. 6(a)] shows that Aao increases from 
close to for W much below TK/bo (indicated by vertical 
dashed lines) to a maximal value of 

^lo = Vi^y + (ly = V2 (56) 

for W ^ Tx/b^. This maximal value reflects the dis- 
placed charge Ach [cf. Eq. (7)] induced by a very strong 
quench: both and are ~ (or ~ 1) if the 
level position is far above (or below) the Fermi energy, 
4 ^ -U/2 + W/2 > (or < 0), cf. Fig. 7(b), thus the 
displaced charge associated with both til and nn is 1. 
(The contribution to Ach from the leads turns out to be 
negligible here,^^ since for sufficiently large W the Fermi 
sea is essentially decoupled from the dot.) 

D. Summary for PS without sensor 

The results of this section can be summarized as fol- 
lows: (i) The energy scale setting the width of PS is pro- 
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portional to Tk- (ii) This can directly be attributed to 
AO: as shown in Fig. 9, the ground states of two configu- 
rations on opposite sides of the switching points exhibit 
strong AO when their level positions differ by more than 
Tk/^o- Thus, quantum fluctuations between them, in- 
duced by operators such as Y and , are strongly sup- 
pressed, (iii) For the present model PS will always be 
continuous as a function of £(j, because (for given U) Tk 
is nonzero for any fixed choice of Fl and Fr (although 
exponentially small), and AO ceases to be important 
(A AO — 0) once comes within Tk/6o of the switch- 
ing point. Conversely, however, it should now also be 
plausible that an essentially abrupt PS will be achiev- 
able if, by a suitable modification of the model, the de- 
gree of AO between the configurations (0,1) and (1,0) 
can be enhanced sufficiently to push Tk to zero even for 
finite Fl and Fr. As pointed out by GBG,^'''^* this can 
be achieved by adding a charge sensor, to which we turn 
next. 
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V. POPULATION SWITCHING WITH SENSOR 



{Ui- Us)/U 




Figure 10: Cartoon of the Hamiltonian (57), describing an 
asymmetric SIAM with an additional sensor lead coupled elec- 
trostatically to the left dot. 

In this section we study the effects of adding an electro- 
statically coupled charge sensor to the model of the pre- 
vious section, as proposed by GBG,"'^'''^^ and analyse how 
this enhances the effects of AO. In particular, we show 
that by increasing the sensor coupling strength (Us), the 
effective Kondo temperature (T^) can be driven to zero, 
implying that population switching becomes abrupt. (A 
study of how additional leads increase the effects of AO 
for static quantities has recently been performed in sim- 
ilar context, involving a multi-lead IRLM.''^ ) 



A. Width of switching regime 

GBG proposed to extend the asymmetric SIAM stud- 
ied above by introducing a third lead as "charge sensor" 
for the left dot (see Fig. 10). For simphcity, it is taken 
to have the same density of states as the other two leads, 
but in contrast to the latter, it couples to the left dot 
only electrostatically (not by tunnelling), with interac- 



Figure 11: (Color online) Population switching for the charge 
sensor model of Eq. (57). (a) nnied) (solid lines) and nL(eii) 
(dashed lines) for several values of Us/U, plotted versus 
(sd + U/2)/Wps in a pseudo-logarithmic fashion ("pseudo" 
in that the x-axis is plotted logarithmic with positive and 
negative values to the left and right of the switching point, 
respectively, represented by the vertical solid line). The hor- 
izontal light solid lines indicate the values of which define 
the widths M^ps of the PS regimes. The noisy behaviour of 
the curves for Us = 5U at small values of Ed indicates that 
our analysis cannot resolve smaller values for Sd as we are 
reaching the limits of double precision numerical accuracy, 
(b) Inset: T^/Tk as function of Us/U, showing the rapid de- 
crease of the Kondo temperature with increasing coupling. 
Main panel: In T| / In Tk versus [Us — Us)/U, plotted on a 
log- log scale (dashed line with dots), together with a linear 
fit using Eq. (59) (solid line). 

tion strength Us (with cs = ^^s): 

H = ^siAM (57) 

e 

A plot of ?iL and riR as functions of Sd for this model 
looks essentially similar to Fig. 7(b), showing population 
switching at Sd = —U/2. However, when the strength of 
the coupling Us is increased, the width of the PS, say 
Wpg, is strongly reduced below the value Wps it had 
for Us = 0, as predicted by GBG. This is illustrated in 
Fig. 11(a), which shows tir (solid lines) and riL (dashed 
lines) as functions of [sd + U /I) /W-ps, using a logarithmic 
scale to zoom in on the immediate vicinity of the PS. In 
fact, as Us approaches a critical value ?7g, the width VFpg 
drops exponentially towards zero, until it becomes too 
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small to be resolved within double precision numerical 
accuracy. 

The behaviour of Wp^ is mimicked by that of the 
Kondo temperature, calculated via the pseudospin sus- 
ceptibility, = 1/xs- We find that it decreases relative 
to its ;7s = value Tk, precisely in proportion to Wpg, 
such that 



2| 
Tk 



(58) 



holds within our numerical accuracy. 

The transition from a continuous to an abrupt PS as 
Us crosses C/g has been predicted to be of the Kosterlitz- 
Thouless type.^^'^* This implies that is expected to 
approach zero according to 



(59) 



where 7 = ^1/2- To test whether our data is con- 
form to this expectation. Fig. 11(b) shows In(T^)/ In(TK) 
vs. (J7g — U) on a log-log plot. Indeed, we find a 
straight line for Us not too close to J7g, consistent with 
Eq. (59). We extract the values 7 = —0.54 ± 0.06 and 
Ug/U = 6.78 ± 0.32, by making linear fits over several 
somewhat different fitting ranges and taking the average 
and standard deviation of the fit parameters as final fit- 
ting results. The relatively large errors of about 10% are 
a consequence of the fact that it is not possible to obtain 
data for Us closer to Ug, since this would drive r| below 
the level of numerical noise. 

We note that analytical calculations based on Refs. 17 
and 18 [using the more accurate criterion, J^iJJ^) ~ 
Jxy[Us) in the notation of these papers] predict the crit- 
ical interaction to be Ug/U ~ 7.6. The agreement of 
this prediction with the numerical result of 6.8 is quite 
respectable, given the inaccuracies in both the numerical 
and analytical calculations [for the latter, inaccuracies 
arise since the cutoff scheme employed in the analyti- 
cal calculation is different from the one realized numeri- 
cally. The cutoff appears explicitly in the arguments of 
the functions Q in Eqs. (6) to (10) of Ref. 17]. 

Though the above results unambiguously show that 
the width of PS decreases exponentially as Us approaches 
a critical value [/g, an analysis based purely on Wps can 
not access the critical point itself or the regime beyond. 
We therefore proceed now with a numerical calculation of 
the dynamics of the pseudospin-flip operator, for which 
we are not constrained to Us < Us- 
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Figure 12: (Color online) (a) The pseudospin-flip spectral 
function A'p'{uj) [cf. Eq. (53)] for the PS model with charge 
sensor, for several values of Us/U, calculated a,t Ed = —U/2: 
plotting A'^ (uj) / A'^ {T§^) versus io/T§^ yields a scaling col- 
lapse. The general shape of the curves is similar to those 
shown in Fig. 8: for u < we find Fermi liquid behaviour, 
~ ui'^ (dotted line), while for lo > Tk each curve shows a 
nontrivial AO power- law, ~ uj-^+'^'^Y [cf. Eq. (54)], exempli- 
fied by the dashed and dash-dotted lines for Us/U = and 
5, respectively, (b) Comparison of the values for r^y expected 
from Eq. (60) (solid line), or extracted from the spectral func- 
tion A'^{io) in the intermediate-frequency regime between Tk 
and the high frequency maximum (crosses). The relative er- 
rors are below 5%, where the errors decrease with increasing 
?7s for similar reasons as in Fig. 8. The light horizontal line 
indicates 77!- = 1. (We were unable to obtain reliable data 
for Us around 7U, presumably because this is too close to 
US-) (c) A^^{uj)/A^^{Tk) versus cj/Tk for Us = 8U. The AO 
power-law behaviour uj~^^^^y extends down to the smallest 
frequencies accessible, illustrating that the crossover scale 
has become undetectably small. 



B. AO in dynamics of pseudospin-flip operator 

The reason for the C/g-dependence of M^ps and is 
that the introduction of the sensor {Us 7^ 0) increases the 
influence of AO in the leads. As pointed out by GBG, 
the scaling dimension of Y acquires an extra contribution 
^Ag due to the sensor lead: 

4 = i(AL + I f + K-Ar - 1)' + h^l , (60) 



where Ag is given by Eq. (38), with C/g replacing U. 
By increasing Us and thereby Ag, it is thus possible to 
drive riy beyond 1. This will render the pseudospin-flip 
operators Y and Y^ irrelevant, thus suppressing quantum 
fluctuations between the (0,1) and (1,0) conflgurations, 
and, hence, pushing down to zero. 

To check this scenario explicitly, we have studied the 
L/g-dependence of 77^ by extracting it from the spectral 
function (w), calculated at the particle- hole symmet- 
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Figure 13: (Color online) AO for the PS model with charge 
sensing. The exponent Aao [extracted from Eq. (1)] (solid 
lines with dots) is shown as function of quench size W jU 
[Eq. (55)], for several values of Us/U , with fixed values of Fr 
and Fl. We see that T^/bo (indicated by dashed vertical lines) 
is pushed to zero as Us increases past Ug/U ~ 6.78. Already 
for Us/U > 6 the curves are essentially indistinguishable, in 
that they do not deviate from their constant value for all 
W^/l7-values accessible to our analysis. For W 3> T§^/bo the 
exponent Aao saturates to a maximal value given by Eq. (61). 
The corresponding values of Ach [from Eq. (7)] are not shown 
but agree with Aao with relative errors of a few percent. 

ric point for several values of Us- The general shape of 
A'p, shown in Fig. 12(a), is similar to that of Fig. 8(a) 
for Us — 0: For frequencies well below T^, A'^{uj) scales 
as w"^, while in the regime of intermediate frequencies, 
^ ^ '^^hc (cf. Sec. IV B), the spectrum shows AO 
power-law behaviour, ^ 0;^^+^^''*'. Indeed, Fig. 12(b) 
shows that the values for ryy extracted from the spec- 
tra (crosses) agree fairly well with those expected from 
Eq. (60). Moreover, for sufficiently large Us/U, the expo- 
nents tjy increase past 1, confirming that the pseudospin- 
flip operators become irrelevant. 



C. AO induced by quench of level positions 

To further highlight the effect of AO on T^, let us con- 
sider again the quench of level position [Eq. (55)] studied 
in Sec. IV C, and repeat the analysis presented there, but 
now for several different values of Us/U. Figure 13 shows 
the results for the exponent Aao • For large values of W 
the AO factor reaches its maximal value 

^aE" = '\Jw' + + . (61) 

This is similar to Eq. (56) for the model without sen- 
sor, but includes the additional contribution Ag [given 
by Eq. (38), with Us replacing U] from the displaced 
charge induced in the sensor lead by the change in local 
occupancy of the left dot from t^l = to 1. 

The most important feature of Fig. 13 is the fact that 
the crossover scale T^/bo (indicated by vertical dashed 
lines) is rapidly pushed to extremely small values as 
Us/U is increased. Indeed, for Us = 8U, which lies be- 



yond the critical value of U^/U ~ 6.78 discussed above, 
Aao is essentially pinned to its maximal value down to 
the smallest values of quench range W that we can access 
numerically. This is consistent with the fact that the cor- 
responding spectral function A'^{lu) at Us — SU, shown 
in Fig. 12(c), shows nontrivial AO power laws down to the 
lowest frequencies accessible, with no trace of a Fermi- 
liquid uj^. This demonstrates very clearly, if somewhat 
indirectly, that the PS will be abrupt for Us > Ug. 

D. Summary for PS with sensor 

Let us summarize the results of this section, by way of 
extending the list of salient points collected in Sec. IV D. 
(iv) The presence of a charge sensor reduces the crossover 
scale T^, which reaches zero at a critical coupling Ug 
[Fig. 11]. (v) This reduction is due to the increased 
effect of AO in the leads, which increases the scaling 
dimension rjy [Fig. 12]; when the latter passes 1 (cor- 
responding to Us — Ug), the pseudospin-flip operators 
become irrelevant and equals zero, rendering the PS 
abrupt, (vi) Correspondingly, for Us > Ug, the spec- 
trum A'^{lj) shows nontrivial AO power-law behaviour, 
^-i+2riY ^ T^g^y (jown to the smallest frequencies ac- 

cessible [Fig. 12(c)], and a low-frequency regime showing 
Fermi-liquid exponents does not exist. 

VI. CONCLUDING REMARKS 

The goal of this paper was to elucidate the role of 
the Anderson orthogonality catastrophe in giving rise to 
anomalous scaling dimensions in dynamical correlation 
functions for quantum impurity models. To this end, we 
have studied several setups involving (interacting) quan- 
tum dots and (non-interacting) leads. The quantum dots 
and leads may be interconnected electrostatically, or also 
through tunnel-coupling. In our analysis we focussed on 
the asymptotic behaviour of various correlation functions 
G{t) and the corresponding spectral functions A{uj) in 
the limit of long times or low frequencies, respectively. 
Their asymptotic behaviour could be understood via a 
generalized version of Hopfield's rule, whose validity was 
checked and confirmed through an extensive NRG analy- 
sis. As a particular application, we performed a detailed 
study of population switching, both without and with a 
third lead that acts as a charge sensor. We confirmed 
a previous prediction^^^^^ that when the charge sensor 
is sufficiently strongly coupled, population switching can 
turn into an abrupt quantum phase transition. 

Aside from presenting a systematic discussion of the 
generalized Hopfield rule, which, hopefully, will be use- 
ful for practitioners in the fields, several general features 
have emerged from our analysis: 

(1) In the context of a local quantum quench of type 
1, where a change of parameters switches the Hamilto- 
nian from Hi to H{, each lead-dot electrostatic coupling 
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Figure 14: Schematic depiction of an equilibrium spectral 
function A'^{u}) for the cases that the local charge relaxation 
rate oj* is (a) larger or (b) smaller than the lead level spacing 
SE. 



gives rise to an AO factor in the ground state overlap 
KCJilGf)!, reflecting a change in the many-body configu- 
ration of the lead when the charging state of the dot is 

modified. This AO factor scales as Nf_i ^ where iV^ 
is the number of electrons in lead /it and the change 
in the scattering phase, divided by tt, in that lead. (AO 
factors from leads that are not interconnected by tun- 
nelling, so that the total charge within each channel is 
conserved, are multiplicative [Eq. (27)].^^ ) 

(2) AO also arises for a type 2 quench, induced by an 
operator that connects initial and final ground states 
|Gi) and \Gf) lying in dynamically disconnected sectors 
of Hilbert space. In particular, AO influences the corre- 
sponding quench spectral function Ax{i^) which scales as 
Axioo) ~ uj-^+^i [Eq. (23)]. For a Hamiltonian without 
tunnelling terms such as the LCM of Eq. (2), the spectral 
function for — thus scales as Ad{i^) ^ uj^^^^''. 

(3) When a type 2 quench has the form of a tunnelling 
operator, = cd\ the asymptotic power law is modified 
to become A^c^ ^ m-i+(Ad+i)" jgq^ (35c)], implying a 
scaling dimension 77^^^^ = 5(^ti + 1)^- For a particle- 
hole symmetric interaction term [as in Eq. (36)], we have 
-1 < Ad < [Eq. (38)], implying that < ry^^t < 1/2, 
thus tunnelling between a dot and a single lead is always 
a relevant perturbation. 

(4) The scaling exponent can be increased, and AO 
strengthened, by coupling the dot(s) to further leads. In 
particular, leads that couple to the dot only electrostat- 
ically (not via tunnelling) contribute AO exponents of 
the form 1 A^ , and thus enhance AO more strongly than 
leads that are tunnel-coupled [cf. point (3)]. In this way, 
the scaling dimension of the tunnelling operator can be 
increased past 1 [cf. Eq. (60)], and tunnelling rendered 
irrelevant. In such a situation, population switching be- 
comes a quantum phase transition, tuned by gate voltage. 

(5) A particularly revealing way of demonstrating the 



effect of AO for population switching is to calculate the 
exponent Aao for a type 1 quench in which the level 
position is abruptly changed from lying above to below 
the PS point (see Figs. 9 and 13, which are analogous to 
Fig. 6(a) for the IRLM). 

(6) In the presence of tunnelling terms of the form 
{c''d + d^c), operators such as = d} , cd^ and c'd^ do 
not induce a quench, since they do not cause a switch be- 
tween disconnected sectors of Hilbert space. Thus, when 
such an operator acts on the ground state, the result- 
ing state will relax back to the ground state over long 
time scales, say t ^ l/i^*, where uj* represents the local 
charge relaxation rate. 

(7) The corresponding equilibrium spectral function 
AP{uj) thus typically shows trivial Fermi-liquid expo- 
nents [e.g. Eq. (39)] in the regime of very small frequen- 
cies, SE < u! u* , where SE represents an infrared 
cutoff such as the level spacing in the lead. (Throughout 
this paper we took 6E ~ 0, since in NRG calculations SE 
can be made arbitrarily small by using sufficiently long 
Wilson chains.) 

(8) In an intermediate frequency regime uj > uj* , the 
equilibrium spectral function A'^{uj) may nevertheless 
contain signatures of anomalous AO exponents, scaling as 
i^-i+2riY Y<]q. (41)], where 77^ represents the scaling 
dimension of Y calculated in the absence of tunnelling. 
Thus, such exponents may be extracted by focussing on 
this regime of intermediate frequencies [as done in Figs. 4, 
8 and 12]. This is schematically indicated in Fig. 14(a). 

(9) If AO can be made so strong that the scaling di- 
mension r/y of the operator is larger than 1, the scale 
a;* is pushed below SE (or, in the context of NRG, below 
the level of numerical noise). In this case, the regime of 
anomalous AO scaling uj~^~^'^^y will extend all the way 
down to the smallest frequencies accessible [e.g. 12(c)], 
as schematically indicated in Fig. 14(b). 

To conclude, we note that cases where AO dominates 
in the low frequency limit such that uj* ~ 0, [as in point 
(9)], quantum fiuctuations of the charge on the dot(s) 
are essentially completely frozen out. At zero tempera- 
ture and in the absence of any extraneous decay mech- 
anism, the system will remain localized in a particular 
local charge configuration. Thus, varying the gate volt- 
age in such a situation may lead to hysteretic behaviour. 
It would be very interesting to experimentally search for 
such signatures of the freezing out of charge fiuctuations 
by performing linear response measurements at the PS 
point. 
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Appendix: Fermi-liquid spectral functions 

In this appendix we study analytically the low energy 
(w < w*) behaviour of the spectral functions of the IRLM 
(Sec. IIIB) and the PS setup (Sections IV B and VB). 

Let us start from the noninteracting resonant level 
(Eq. (36) with U = 0). In that case an elementary cal- 
culation gives for the retarded dot Green function,'*^ 



^-^d- :^,Gc ■ (w) 



1 



(A.l) 



where Q^'^ is the retarded c Green function for F = 0, 
and we assumed the wide band limit (used just to sim- 
plify expressions, but actually not essential for any of the 
following arguments) Q^''^{ijj) = —iirp = — i7r/2 in units 
where D = 1. The imaginary part of the retarded Green 
function gives (up to a factor of — I/tt) the well-known 
Lorentzian spectral function 



A7{u) = 



1 



r 



(A.2) 



Thus, at low energies (w <C F) A'^^{ijo) becomes a con- 
stant, corresponding to nfj^^ = 1/2 [which reproduces 
Eq. (39a)]. This behaviour is easy to understand: In the 
absence of tunnelling ^°'^'°(a;) = p is constant, reflecting 
the constant local density of states of the lead electrons 
near the end of the lead. In the presence of tunnelling, 
at low energy the dot level is well-hybridized with the 
lead, and assumes the role of the end of the lead, thus 
featuring the slowly-varying low-energy spectral function 

Based on similar arguments, one would expect that, 
in the presence of tunnelling, A'^^{lli) is still constant at 
low-energies, since in that limit the small spatial separa- 
tion between the dot and the end of the lead should be 
unimportant. However, commensurability at half filling 
(particle-hole symmetry) makes things bit more compli- 
cated. An explicit calculation gives: 



-nrp 



^ - grf 
uj — ed + iT 



(A.3) 



Thus, when Ed is nonzero, we indeed get a constant low 
energy limit, i.e. 77°'^ = 1/2. However, when Ed — 



(the value used throughout this paper for the IRLM), 
g^{uj) ^ w while A°^{uj) ~ w^, corresponding to ry^'^ = 
3/2. To understand this behaviour, let us examine a 
half infinite tight-binding chain with lattice spacing a and 

Hamihonian 7?tb = Er=i(*Ui*n + h.c). Taking the 
continuum limit in the standard way, we can expand the 
fast-varying annihilation operators in terms of slowly- 
varying (on the scale of the Fermi wavelength) right /left 
moving fields tpR/Lix), with x ~ na: 



. ik-pna 



ipR^na) 



ipL{na) , (A.4) 



where fcp is the Fermi wavevector. From the boundary 
condition ^fg = one gets "01.(0) = — V'fl(O): so we can de- 
fine the single slowly-varying field i^i^) by ipi^) = '4'r{x) 
if x > and ^(x) = —^l{—x) if a; < 0. Then: 



V.(- 



(A.5) 



At half filling, kpa — tt/2, we get at the site next to the 
boundary 

*„=2 - -^(2a) + 0(-2a) - -4a9,V(0) , (A.6) 

The same thing happens at the first site (n = 1) when 
we attach a dot, since at low energies the dot behaves 
as the new first site. The spatial derivative is equiva- 
lent to a time derivative, up to the Fermi velocity vp. 
This extra time derivative is responsible for the vanish- 
ing of the spectral function Al'^{u}) for a; — >■ 0. Since 
we have derivative for both c and in the Green func- 
tion, and each gives an extra factor of we end up with 



This behaviour depends on being at half 



filling (particle-hole symmetry), hence is modified when 
Ed is not zero. 

Now we can discuss the higher spectral functions, 
A'7.f{uj), and A'7{uj). These are the imaginary parts of 
the corresponding retarded Green functions, up to a fac- 
tor of — I/tt. The retarded Green functions are in turn the 
analytical continuation of the temperature Green func- 
tions to the real frequency axis. And the temperature 
Green functions can be found in the noninteracting case 
using Wick's theorem.*^ 



Performing these calculations for A^J^^^icd), one gets 



In 



1^ - £d 



iV 



(A.7) 



uj{u + 2ir) 



In 



Ei -iiu + ivy^ 



Concentrating on w <g; F one finds rj^ 



r2 



for Ed ^ and 77^' 



eq 



v. 



eq. 



'S'l + - 1 



2 for Ed — [the data in 



Fig. 4(c) corresponds to the latter case, which reproduces 
Eq. (39c)]. This simple summation of scaling dimensions 
is natural here, since there is only one possible different- 
time Wick-pairing, of each single-particle operator with 
its conjugate. 

For ^^^(oj), however, there are two different-time 
Wick-pairings, causing cancellations, and resulting in: 



P^ 



2Ed 



2Ed 



2iF, w 
In 



£d + iT 



~Sd + iT 



(A., 
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Concentrating again on w <C F one finds now that 77^^ = 2 
for all values of Ed [which reproduces Eq. (39c)]. The 
reason is that in the low-energy continuum limit the 
product dc becomes the product of annihilation opera- 
tors at almost the same point. Hence, one should expand 
in the distance between d and c (lattice spacing). The 
leading term (with no spatial derivatives) vanishes by 
Pauli's principle; the next term involves a spatial deriva- 
tive, leading to a factor of uj, similarly to the arguments 
above. Another factor of uj comes from the operator c^d} 
appearing in the definition of Q^'ci'^)- Thus, at low ener- 
gies we end up with A'^J^{uj) ^ uj^ even for Sd 7^ 0. 

Although the above calculations were performed for 
the noninteracting case, the qualitative arguments ex- 
plaining the low-energy behaviour are valid even when 
U > 0. Moreover, since the system flows to the same 
fixed point for all values of [/ > 0, the low energy power- 
laws are in any case independent of U. Our numerical 
results (Fig. 4) are in agreement with this picture. 



Let us now turn to the low-energy behaviour of the PS 
setup in the case where PS is continuous. At low-energy 
the system is governed by Kondo physics, Eq. (49),^^"''® 
where the L-R degree of freedom plays the role of a 
pseudo-spin. The equivalence to Kondo continues to 
hold even in the presence of a charge sensor, as shown 
by GBG}'''^^ The operator Y = dRd-L^RC^L (which is 
relevant in the continuous-PS phase) is the pseudospin- 
fiip local exchange term between the dot and the lead. 
Similarly, the IRLM is also equivalent to the Kondo 
model, ^^'^^''^^ with the role of spin replaced by the charg- 
ing state of the dot. The pseudo-spin local exchange term 
is simply dc' . Hence, when the parameters are properly 
mapped, the spectral functions A'^{lo) and y^°^^(u;) are 
equivalent when the Kondo description applies (i.e. for 
CO <^ D for the IRLM, and uj <^ Whe for the PS setup). 
In particular, AP{uj) should exhibit an behaviour at 
low energy, similarly to A'^^^ (w) for Sd = 0, as the NRG 
data shows [dotted line in Fig. 8(a) and Fig. 12(a)]. 
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